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We ascertain, following ideas of Arnesen, Bose, and Vedral 
concerning thermal entanglement [Phys. Rev. Lett. 87 
(2001) 017901] and using the statistical tool called entropic 
non-triviality [Lamberti, Martin, Plastino, and Rosso, Phys- 
ica A 334 (2004) 119], that there is a one to one corre- 
spondence between (i) the mixing coefficient a; of a Werner 
state, on the one hand, and (ii) the temperature T of the 
one-dimensional Ifeisenberg two-spin chain with a magnetic 
field B along the axis, on the other one. This is true for 
each value of B below a certain critical value Ba. The perti- 
nent mapping depends on the particular B— value one selects 
within such a range. 

Pacs: 03.67.-a; 89.70.-)-c; 03.65.-w; 02.50.-r 



I. INTRODUCTION 

Entanglement is one of the most fundamental issues of 
quantum theory [1-6] and the so-called Werner states [7] 
have played a distinguished role in the unravelling of the 
fascinating issues at play (see, for instance, [8-10]). The 
Werner density matrix reads 



pw = 



1 



(1) 



where [$+) is a Bell state (maximally entangled). The 
state (1) is separable (unentangled) for the mixing coef- 
ficient a; < 1/3 [7]. For x > 1/3 they are entangled and 
violate the CHSH inequality for x > 1/^2 [8-10]. We see 
that Werner states are mixtures of noise and a maximally 
entangled state, and therefore, for values of the mixing 
parameter x > 1/3 they are entangled, violate the CHSH 
inequality and exhibit nonclassical features [8-10]. 
Following the interesting work of Arnesen et al. [11], we 
concern ourselves with the issue of thermal entanglement 
and consider the Hamiltonian for the ID Heisenberg spin 
chain with a magnetic field of intensity B along the z-axis 



N 



(2) 



where al. y ^ stand for the Pauli matrices associated with 
spin i and periodic boundary conditions are imposed 



AT -1-1 



^u)- Jh gives the strength of the spin-spin re- 



pulsive interaction (only the anti-ferromagnetic {Jh > 0) 



instance is discussed). If we limit ourselves to the case 
A'^ = 2, we will be dealing with two spinors, i.e., with a 
two-qubits system. For "thermal equilibrium" we should 
consider [11] the thermal state 



exp(-fc^) 
ZiT) 



(3) 



with Z{T) the partition function. Expressing both H 
and p{T) in the computational basis |00), |01), |10), |11) 
we obtain 



H = 



f 2Jh + 2B 

-2Jh 4:Jh 

AJh -2Jh 
V {) 2Jh -2B ) 



(4) 



After defining, for convenience's sake, 
ewray = exp {~2w - 2y); 
e-wp = exp (— 2w) + exp (6w); 
eyjra = exp (-2w) - exp (6w); 



e-wpy = exp {-2w + 2y), 
with w = Jn/kBT and y 

/ ^wmy 



B/ksT, we also get 
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The concurrence of p(T) reads [11] 



(5) 



C = 0; 

C = 



1 -1- e-'^y + e^y + 



for T > Te, 

forT < Tr 



(6) 



For our purposes we must emphasize that there is no en- 
tanglement beyond a certain critical temperature Tc = 
8JH/(fcBln3) [11]. Clearly, the temperature plays a sort 
of "mixedness-role" , since it is well-known that entangle- 
ment vanishes for a high enough degree of mixing. Notice 
that Tc is independent of B and thus an intrinsic struc- 
tural property that depends (linearly) just on the interac- 
tion strength J^. For strong enough coupling strengths 
we encounter non null entanglement for a wide range of 
temperature values. 

Also, there is a change in the structure of the ground 
state of hamiltonian (2) when the magnetic field reaches 
the critical value Be = 4J/f. 
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Remarkably, the Werner states (1) and the thermal states 
(3) of the ID Heisenberg model for iV = 2 can be re- 
lated by a one to one correspondence between the mix- 
ing parameter x and the temperature T for each value of 
B < Be- By means of this relation we explore for both 
states the evolution of the entanglement of formation as 
a function of the temperature. 

The paper is organized as follows: In Section II we 
introduce the entropic non triviality measure based on 
the Jensen-Shannon divergence, which is the basic tool 
that provides a deeper insight into the behaviour of the 
entanglement of formation as a function of T around Be- 
Our results are reported in Section III, where we provide 
the exact mapping between Werner states and the afore- 
mentioned thermal states. Finally some conclusions are 
drawn in Section IV. 



II. THE JENSEN-SHANNON DIVERGENCE 
AND THE ENTROPIC NON-TRIVIALITY 
MEASURE 

We review now a statistical information measure with 
which we will be concerned in the present work. Let 
P(k) G C TZ^ , with k = 1,2, denote two different prob- 
ability distributions for a particular set of N accessible 
micro-states. The components of the two probability vec- 
tors P(fc) must satisfy the following two constraints: a) 



pf^ = 1, and h) Q< pf < 1 Vj. The set n de- 



fined by these constraints is the simplex Sjv, which is a 
convex (TV — l)-dimensional subset of R,^ . A quite im- 
portant, information-theoretical based divergence mea- 
sure between P(i) and P(2) was originally introduced by 
Rao [12] and used by several authors [13]. It came after- 
wards to be called the Jensen-Shannon divergence (JSD) 
[14-16] that 

• induces a true metric in O c TZ^ , being indeed the 
square of a metric [16] , and 

• is intimately related to the KuUback-Leibler rel- 
ative entropy K for two probability distributions 
P(i) and P(2), given by [17] 

^[^(1)1^(2)] = E^^r )• (7) 



We first define 

^0[P(l),P(2)] = ^[^(1)1(^^(1) + ^4))]. (8) 

and then the symmetric quantity 

MP{1),P{2)] = MPil),Pi2)] + MP{2),Pil)] (9) 
= 25[^P(i) + ^P(2)] 



-5[P(i)]-5[P(2)], 

where S = — '^jPj Inpj is the Shannon logarithmic in- 
formation measure. Let now tti, 7r2 > 0; tti -|- 7r2 = 1 
be the "weights" of, respectively, the probability distri- 
butions P(i), P(2). The JSD reads 



-P(i),-P(2) 7riP(i) +7r2P( 



(2) 
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(10) 



which is a positive-definite quantity that vanishes iff 
P(i) = P(2) almost everywhere [14,15]. In the particu- 
lar case to be used in this work, tti = 772 = 1/2, the 
measure (10) is symmetric. Notice also that J2^2 = 1 Jj^. 

Another statistical tool that we need is the so-called 
entropic non- triviality measure [18]. The statistical char- 
acterization of deterministic sources of apparent random- 
ness performed by many authors during the last decades 
has shed much light into the intricacies of dynamical be- 
havior by describing the unpredictability of dynamical 
systems using such tools as metric entropy, Lyapunov 
exponents, and fractal dimension [19]. It is thus possible 
to i) detect the presence and ii) quantify the degree of 
deterministic chaotic behavior. Ascertaining the degree 
of unpredictability and randomness of a system is not 
automatically tantamount to adequately grasp the corre- 
lational structures that may be present, i.e., to be in a 
position to capture the relationship between the com- 
ponents of the physical system. Certainly, the opposite 
extremes oi i) perfect order and ii) maximal random- 
ness possess no structure to speak of. In between these 
two special instances a wide range of possible degrees of 
physical structure exists, degrees that should be reflected 
in the features of the underlying probability distribution 
(PD). One would like that they be adequately captured 
by some functional of PD in the same fashion that Shan- 
non's information measure captures randomness. A can- 
didate to this effect has come to be called the entropic 
NON-TRIVIALITY (also, "statistical complexity") C [18], 
that should, of course, vanish in the two special extreme 
instances mentioned above. 

In order to illustrate these assertions we can refer to 
the celebrated logistic map P : a;„ — > Xn+i [20], where 
one focuses attention upon the ecologically motivated, 
dissipative system described by the first order difference 
equation 

Xn+i = r Xn (l-Xn) ( < a;„ < 1 , < r < 4 ) (11) 

whose dynamical behavior is controlled by r. For val- 
ues of the control parameter 1 < r < 3 there exists 
only a single steady-state solution. Increasing the con- 
trol parameter past r = 3 forces the system to un- 
dergo a period-doubling bifurcation. Cycles of period 
8, 16, 32, • • • occur and, if r„ denotes the value of r 
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where a 2" cycle first appears, the r„ converge to a Hm- 
iting value Too — 3.57 [20]. As r grows still more, a quite 
rich, and well-known structure becomes apparent. A cas- 
cade of period-doubling occurs as r increases, until, at 
Too, the maps becomes chaotic and the attractor changes 
from a finite to an infinite set of points. For r > r^a 
the orbit-diagram reveals an "strange" mixture of order 
and chaos, with notable windows of periodicity beginning 
near r = 3.83. As shown in [15], the measure C is able de 
give detailed account of the pertinent dynamical features 
(see, for instance. Figs. 2-4 of [15]). 

Typically, the measure C is the product of two quan- 
tities: (1) a normalized entropy (i.e., whose values range 
between and 1) H = S/Smax (the denominator is the 
largest possible value for S) and (2) a chosen "distance" 
d{P, Pu) in probability space that measures "how far" (in 
this space) the actual PD P lies from the uniform PD (of 
maximal entropy) P„, i.e.. 



C==Hd{P,Pu) 



(12) 



(see, for example, [18,21 26] and references therein). In 
the present work we are concerned with the probability 
set (density matrix' eigenvalues) associated to any state 
p and then, with reference to (12), take as H the nor- 
malized von Neumann entropy Hy^ of the relevant state 
p and as d the Jensen-Shannon divergence that will here 
measure the "distance" from the relevant state to the 
maximally mixed one pmm (proportional to the identity 
matrix I). Our entropic non- triviality acquires then the 
aspect (Cf. Eq. (10)) 



Cjsip) = J^'^'^'\p,Pmm)H,n{p). 



(13) 



One may wonder whether it is possible to built up an en- 
tropic non-triviality measure that would depend not on 
density matrices but directly on wave functions. The con- 
comitant extension is of a rather non trivial nature, but 
can be accomplished and steps toward its implementation 
are currently being undertaken. They will be published 
in the near future. 



III. RESULTS 

Our main interest refers to the "temperature- 
evolution" of both Cjs and the entanglement of forma- 
tion Ef as T diminishes, and, eventually, in the limit case 
T — > 0, where a change in the structure of the ground 
state of hamiltonian (2) is detected for = ^Jh- For 
convenience, wc take Jh = 1 from now on. In the limit 
T — » one can state that 

• \) [B < T = {)] the gromid state (gs) is non- 
degenerate and equal to the singlet state {Ef — 1, 
Cjs = for all B) 



p{\T 



0) 
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2) [B = Be, T = Q] the gs state becomes two- 
fold degenerate, the corresponding two eigenstates 
being the singlet state and |11) {Ef = Ef, Cjs = 



pf{T = Q) 
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• 3) [S > Be, T = 0] the ground state is now jll), 
which has no entanglement {Ef = 0, Cjs = for 
all B) 



pr(r = o) = |ii)(ii| 



/O \ 





\ 1 y 



(16) 



It is clear that, for fixed Jh and at T = 0, the entropic 

no-triviality is zero for all B except for Be, which is tan- 
tamount to asserting that Cjs "detects" the "critical" 
magnetic strength at which the gs-structure changes. 



We return attention now to the Werner state 



Pw 



|$+)($+| + {l-x)p 



MM, 



(17) 



where 1$"*") is a (maximally entangled) Bell state. This 
state is unentangled for a; < | [7]. Fig. 1 depicts Cjs{p) 
vs. T and Ej{p) vs. T for the ID Heisenberg model in 
the case B = Q. For fixed values of Jr and B, there 

r-^ above which the 



8J 



ks In 3 



exists a critical temperature Tc 
system is no longer entangled. The same quantities are 
depicted in the inset for the Werner state pw If one 
considers the plot Cjs{pw) vs. Ef{pw) and its thermal 
Heisenberg counterpart Cjs{p) vs. Ef{p) for B = 0, the 
two graphs coincide. In point of fact, there is a mapping 
between both states, if we take 



X 



1 -I- e-'^y + e^y + e 



8w 



1 



(18) 



there is a one to one correspondence between the mixed 
parameter x of the Werner states and the temperature T 
of the states (3) for each value oi B < Be (see Fig. 2) so 
that x can be regarded as an effective temperature Teff. 
We remark that 

1. at X = 1 (T = 0) both states are maximally en- 
tangled (and have zero complexity) until a critical 
point is reached, 
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2. at the critical point x = ^ Tc = 3 the 
two states have the same Cjs- Notice that all the 
curves of Fig. 2 (for any several B— values) inter- 
sect at the point {x = 1/3, l/fcsTc), 

3. when X = (T ^ 00) both states are unentangled 
(also with Cjs = 0). 

In the vicinity of Be, the plots Cjs vs. Ef (as pa- 
rameterized by T) provide a good insight into the en- 
suing transition mechanism. In Fig. 3 we plot Cjs{p) 
vs. Ef{p) for several B— values near the critical point 
Be = 4. Note that, as we approach the critical point 
from either the left (inset) or the right, the cntropic 
non-triviality augments. Pay special attention to the 
rightwards (Bf) approach. We not only appreciate a 
maximum in Cjs, but notice also that Ef is optimal at 
B = 4.1. Thus, maximum Cjs equals maximum Ef at 
the critical point. Vertical and horizontal dashed lines 
correspond to critical values at T = 0,B = Be- As 
seen from p^*, there is sudden entanglement-change as 
B crosses B = B^- As Cjs can only detect changes in 
state-structure, its variations can only be produced gs- 
transitions. Its "transition-detection" feature improves 
its accuracy in the limit T ^ 0. 

Transition details near Be = A are examined in an even 
closer fashion in Fig. 4, a Cjs vs. Ef plot. We depict 
things for B+ = 4.001 (soUd line) and for B' = 3.99 
(dot-dashed line). Zone I is the region with Ef > Ej^ 
{B~ ) , wherefrom, as we increase T, the state loses entan- 
glement while increasing its entropic non-triviality (dot- 
dashed line) . See in the inset (solid thick line) the cross- 
ing at Ej , enhanced in order to appreciate details of the 
evolution towards maximal entanglement and maximal 
complexity in the limit T ^ 0. On the other hand, if the 
initial point occurs at B+, no matter whether the tem- 
perature changes start at either T = or T = c», we end 
up at the same Ef = point. In particular, for B = Be, 
there is a single curve that connects the T = critical 
point with the T — > c» limit (where Cjs = Ef = 0). 
The double-arrow connecting the curves in the two zones 
illustrates the fact that a path beginning in Zone I (dot- 
dashed line) cannot smoothly merge with the curve de- 
picted in Zone II (solid line), appearances notwithstand- 
ing. A glimpse at the inset clarifies the situation. Vertical 
and horizontal dashed lines correspond to critical values 
at T = 0, B = Be- 

Fig. 5 is a Cjs vs. Ef plot that intends to de- 
pict things in the immediate vicinity of the critical point 
Be = 4. As in Fig. 4, here we fix a very low T (of the 
order 10~^) while varying B. The solid line corresponds 
to the entropic non-triviality. See the quasi-circular "mo- 
tion" around Be as we decrease the temperature. Inset a) 
depicts Cjs vs. Ef, which, surprisingly, coincides with 
the "envelope-curve" (Zone I -|- Zone II) of Fig. 4. This 
remarkable fact can be explained because in Fig. 4 both 



the "upper" branch of the curves in Zone I and the curve 
of Zone II are drawn for low T-values nearby B = Be. 
Since here we "sit" at the critical point (and very low 
temperatures arc used) both curves of inset a) look like a 
continuous lines. In other words, only a low S— range of 
values around Be is relevant at such low temperatures. 
Inset b) depicts the behavior of the degree of mixture 
R = l/Tr[p{Tf] vs. Ef. In the limit T 0, R reaches 
R = 2, as expected. Vertical and horizontal dashed lines 
correspond to critical values a,t T = 0,B = Be- 

IV. CONCLUSIONS 

In this paper we have shown that the Jensen-Shannon 

cntropic non-triviality measure Cjs is able to detect the 
"critical" magnetic strength at which the gs-structure 
changes. Also, via Cjs we discover that there is a one 
to one correspondence between Werner states and the 
Heisenberg anti-ferromagnet for values oi B < Be- As 
shown in Fig. 2, it is possible to assign to the mix- 
ing Werner parameter x an effective temperature Teff 
for each value oi B < Be- Finally, the physics in the 
immediate vicinity of the critical point Be = 4 can be 
analyzed in detailed fashion using Cjs- In the one di- 
mensional Heisenberg model, for two spinors, maximum 
entanglement implies maximum entropic non-triviality in 
the vicinity of the corresponding critical point. 

Acknowledgements 

This work was partially supported by the MEG grant 
BFM2002-03241 (Spain) and FEDER (EU), by the Gov- 
ernment of Balearic Islands and by CONICET (Argen- 
tine Agency). 

FIGURE GAPTIONS 

Fig. 1- Evolution of (1) the entanglement of formation 
Ef and (2) the Jensen-Shannon complexity Cjs vs. the 
temperature T for the thermal state p{T) of two-qubits 
in the ID Heisenberg model. For given values of Jh 
and B there exists a critical temperature Te = ^.^'"^jgg 
above which the system is no longer entangled. The same 
quantities are depicted in the inset for the Werner state 
pw- There is a one to one mapping between both type 
of states for B < Be. See text = for details. 

Fig. 2- Mapping of the mixed Werner parameter a; as a 
function of 1/T for several values of B. The horizontal 
lines correspond to a; = 1/3 and x — 2/3, respectively. 

Fig. 3- Plot of the complexity Cjs vs. the entanglement 
of formation Ef near the critical point Be = 4. As we ap- 
proach the critical point from either left (inset) or right, 
the complexity augments. In the case we not only 
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have a maximal Cjs, but also an optimal Ej. Vertical 
and horizontal dashed lines correspond to critical values 
at T = 0, B = Be- See text for details. 



Fig. 4- Same as before, but with additional details. We 
have B+ = 4.001 (solid line) and B' = 3.99 (dot-dashed 
line). Zone I (II) denotes the region with Ef > and 
B~ {Ef < Ej and -B+). In the inset the crossing re- 
gion at Ej is enhanced in order to see how details of the 
evolution around maximum cntangiomont and maximum 
complexity in the limit T ^ 0. See text for details. 

Fig. 5- Plot of the complexity Cjs in the vicinity of the 

critical point Be = 4. The inset a) depicts Cjs vs. Ef, 
which, surprisingly, coincides with the enveloping curve 
(Zone I + Zone II) of Fig. 3. Inset b) depicts the be- 
haviour of the degree of mixture R = l/Tr[p(T)'^] vs. 
Ef. In the limit T — > 0, i? reaches the value R = 2, 
as expected. Vertical and horizontal dashed lines corre- 
spond to the critical values at T = 0,B = Be- See text 
for details. 
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